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ABSTRACT 

OH is an important molecule in the H2O chemistry and the cooling budget of star-forming regions. The goal of the Herschel key 
program 'Water in Star-forming regions with Herschel' (WISH) is to study H2O and related species during protostellar evolution. Our 
aim in this letter is to assess the origin of the OH emission from star-foming regions and constrain the properties of the emitting gas. 
High-resolution observations of the OH J = 3/2 - 1/2 triplet at 1837.8 GHz (163.1 pm) towards the high-mass star-forming 
region W3 IRS 5 with the Heterodyne Instrument for the Far-Infrared (HIFI) on Herschel reveal the first hyperfine velocity-resolved 
OH far-infrared spectrum of a star-forming region. The line profile of the OH emission shows two components: a narrow component 
(FWHM w 4 - 5 km s"') with partially resolved hyperfine structure resides on top of a broad (FWHM a; 30 km s"') component. 
The narrow emission agrees well with results from radiative transfer calculations of a spherical envelope model for W3 IRS 5 with a 
constant OH abundance of xqh ~ 8 x 10"'. Comparison with H2O yields OH/H2O abundance ratios of around 10"' for T > 100 K and 
around unity for T < 100 K, consistent with the current picture of the dense cloud chemistiy with freeze-out and photodesorption. The 
broad component is attributed to outflow emission. An abundance ratio of OH/H2O > 0.028 in the outflow is derived from comparison 
with results of water line modeling. This ratio can be explained by a fast J-type shock or a slower UV-irradiated C-type shock. 

Key words. Astrochemistry — Stars: formation — ISM: molecules — ISM: jets and outflows — ISM: individual objects: W3 IRS 5 
1 . Introduction iGoicoechea et al.ll2006l) and with PACS ('Wamofler et al."201 



/-^ XT , • ■ on board the i/ersc/ie/ Space Observatory (Pilbratt et al. 201 „j,. 

^ Newly-formed stars inject large amounts of energy into the am- previously, OH was mainly observed through maser lines at 

^ bient interstellar material through shocks and radiation, heat- cm wavelengths (.Gaume & Mutill [19871 for models see e.g. 

^ ing the suiTounding gas and dust In these warm regions, H2O |Cesaroni & Walmslevlll991h . Observations carried out with ISO 
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becomes one of the most abundant gas-phase molecules be- p^^^ ^^^^ ^ ^^^^^1 resolution needed to distinguish 

cause water ice evaporates from the grain mantles and gas-ph^^^^^ ^^^^^^^ Q^^ ^^^..^^^ the quiescent envelope and the 

formation routes become available^ The hydroxyl radical (OH) ^^^j^^^^ ^.^^^ information on the line shape is obtained. 

IS closely hnked to both the H2O formation and destruction p^rthermore, the hyperfine structure of the OH lines remaines 

through OH + H2 « H2O + H as well as a byproduct of the H2O unresol ved. The Heterodyne Instrument for the Far-Infrared 

photodissociation process m the presence of UV photons. A sig- ^j^jpj^ |de Graauw et alJ iM on Herschel offers an adequate 

nificant fraction of the gas cooling occurs through line emission ^ ^^^^^ resolution to resolve the 1837.8 GHz triplet of OH. 

in the far-mfrai-ed including lines of [O i], [C n , CO, H2O, and ^j^^^^ j^^j ^ ^-^^ ^-^^^^ constrain the tempera- 

OH (Giannim et al. 2001; van Kempen et al. 2010). ^^^^ ^^^-^^^ ^ ^j^^^^^^ 

Observations of the OH far-infrai-ed transitions first be- qh triplet was previ ously detected with HIFI from Orion KL 

came possible w i th the Kuiper Airborne Observatory (e.g . Jr'r^^Uott at 0I lhni nK 

IStorev et all [19811 Ifietz & Boreikd [19891 iMeMck et al.lll990[ ICrockett et al. || 2010D . 
and references therein). OH was fr equently d etected in 

star-forming regions with ISO (e.g. iGiannini et al. 200il This letter presents the first hyperfine velocity resolved ob- 

servations of OH obtained using HIFI from the well-studied 

* Herschel is an ESA space observatory with science instruments high-mass star-forming region W3 IRS 5 located at a distance 

provided by European-led Principal Investigator consortia and with im- of 2 kpc ((Hachisuka et al. 200 61) with a total luminosity of 

portant participation from NASA. ~ 10^ L0 (e.g. iHelmich & van Dishoecklll997HBoonman et alJ 
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2003 : Ivan der Tak et al l2005t IChavama et al. l2010t iBenz et all Table 1. Molecular data of the observed OH transitions between 
2O10l) . the 3/2 and 1/2 excited states. 



2. Observations and data reduction 

The OH triplet at 1837.747, 1837.817 and 1837.837 GHz from 
W3 IRS 5 (a2000 = 02h25'"40;60, ^2000 = +62°05'51'.'0) was ob- 
served with HIFI on Herschel. The observations were carried out 
as part of the 'Water in Star-Forming Regions with Herschel' 
(WISH) key program (Ivan Dishoeck et al.ll201 1I) . W3 IRS 5 was 
observed on July 29* 2010 (obsid 1342201666) in dual beam 
switch mode with an on-source integration time of 17 min. The 
beam size (HPBW) at 1837 GHz is about 12". The average 
system temperature was 1245 K. The wide band spectrometer 
(WBS) offers a nominal resolution of 1.1 MHz, corresponding 
to a velocity resolution of about 0.18 km s"' at 1837 GHz. The 
OH triplet at 1834.7 GHz was not observed. 

HIFI data were reduced using t he Hersch el interactive pro- 
cessing environment (HIPE v4.0.0, IOtil2010l) and the GILDAS- 
CLASSQ software. The H and V polarizations were combined, 
reaching Tnns » 0.1 K . We subtracted a first order polynomial 
from the spectra and calibrated to r^b scale using a main beam 
efficiency of 0.70. Molecular data (Table [TJ are taken from the 
Leiden atomic and molecular database (Schoier et al. 

12005 ) with data from the JPL catalog] (Picket t et alJll 9981) and 
collisional rate coefficients bv .Off'er et al. (1994k 

The calibration uncertainty of band 7b is currently estimated 
to < 32% and is mainly caused by the unknown side band gain 
ratio. The uncertainty reduces to about 10% for a ratio close to 
unity. The errors on the integrated intensities are only fitting er- 
rors and do not include the calibration uncertainty. The velocity 
calibration is uncertain by ~ 0.2 km s"' and could partially ex- 
plain the ~ 0.5 km s"' offsets of the line peaks from the source 
velocity of Ukr — -38.4 km s"' (.van der Tak et al..,2000i) . 

3. Results 

HIFI clearly detects the OH triplet at around 1837.8 GHz 
(163.1 yum) towards W3 IRS 5 (Fig. [TJ- For the first time, the 
hyperfine components of the line triplet are spectrally resolved. 

Two components dominate the line shape: the line profile 
shows narrow components residing on top of a broad emis- 
sion feature, similar to that found by Chavania et al. (2010) for 
H2O lines in W3 IRS 5 and iKristensen et aL (,2010) in low-mass 
YSOs. The narrow components with full widths at half maxi- 
mum (FWHM) of ~ 4 - 5 km s"' are centered close to the po- 
sitions expected from the hyperfine pattern. The lines are split 
by 20.2 MHz (3.3 km s"') towards the blue and 70.2 MHz 
(11.5 km s ') towards the red relative to the middle line. 
Therefore, the {F ^2+^1-; Tab. [B and (F = 1h- ^ 0-) 
transitions are blended, while the {F - \+ — > 1-) transition 
is resolved. The underlying broad component with a FWHM of 
~ 30 km s ' consists of the three blended hyperfine components, 
each with Ay ~ 20 km s The hyperfine pattern is unresolved 
in the broad component. 

Gaussian fits to the narrow components are presented in the 
online appendix (Sec.|A]i. The derived hyperfine intensity ratios 
(1:3 ± 0.5:1.5 + 0.3) deviate from the prediction of LTE in the 
optically thin limit (~ Auigu, 1:5:2). Pumping by far-infrared 
dust continuum radiation, line overlap, and optical depth may 

' http://www.iram.fr/IRAMFR/GILDAS 
^ http ://w w w. strw. leidenuni v. nl/ ~moldata/ 
^ http://spec.jpl.nasa.gov 
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Fig. 1. HIFI spectrum of the OH triplet at 1837.8 GHz with half 
the dual-sideband continuum. The expected positions of the lines 
and the source velocity are labeled. The blue lines show the best 
fit from the slab models and the outflow component separately. 



all contribute to the observed hyperfine anomaly. Modeling re- 
sults (Sec. 13. Il l indicate that at least the 2-1- —> 1- line is opti- 
cally thick. Based on comparison with spherical envelope mod- 
els (Sec. 13.2b . we attribute the narrow components to OH emis- 
sion from the hot parts of the envelope. 

The broad component in OH is most likely associated with 
outflow emission. Its integrated intensity is around 37 K km s"' 
(Sec. |A|. A similar component can also be identified in the 
line profiles of other species like CO and H2O (Sec. |B]i from 
rece nt Herschel resul ts and ground-based o b servations (e.g. 
Has egawa et al.1 11994 iBoonman et"al] 120031: IChavama etakl 
2OI0I) ! 

3.1. Slab Models 

As a first modeling step to derive OH colu mn densities, the slab 
method outlined in iBruderer et al.l (1201 Oi appendix B) is used. 
This method calculates the molecular spectrum from two slabs 
- one representing the envelope and one the outflow - in front 
of a continuum source, determined from the observed contin- 
uum flux. The slabs are assumed to cover the entire continuum 
source. The free parameters describing each slab are the OH col- 
umn density A^oh. the excitation temperature Tex (assumed to be 
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Table 2. Beam-averaged column density A^oh, excitation tem- 
perature Tex, line width Au, and line position uisr of the best fit 
slab model. 



A'OH ^ 
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^ ex 


Av 




fit range 


[cm"-] 


[K] 


[km/s] 


[km/s] 


[km/s] 


>7.2(13) 


270.0° 


21.9 


-35.7 


[-60,-47.5], [-22,-10] 


E 3.9(16) 


34.5 


3.3 


-37.7 


[-43,-35], [-30,-23] 



" fixed (see text). A{B) = Ax 10®. O = outflow, E = envelope. 



equal for aU hyperfine transitions), the line width Au, and the po- 
sition of the line center t^isr. Overlap between different hyperfine 
components of OH is taken into account. To constrain the free 
parameters (A^oh, ^ex, Av and ui^r), the between observation 
and model is minimized using an uncertainty of T^ms ~ 0.1 K. 

The resolved hyperfine triplet structure of the narrow (enve- 
lope) component allows us to constrain the column density and 
the excitation temperature of the emitting gas simultaneously. 
Results of the best fitting slab models for the narrow compo- 
nent with a reduced = 1-2 are presented in Table |2] A peak 
line opacity of t = 1.8 is reached (2-1 — > 1- transition) and the 
In — > 0- and 1 h — > 1 - peaks have optical depths of r = 0.8 and 
T = 0.4, respectively. Figure ICTI of the online appendix shows 
the Icr and 2cr contours for the column density and excitation 
temperature. The OH column densities within the Icr interval 
range from A^oh = 2 x 10'^ - 1 x 10'^ cm"^, depending on Tgx, 
and thus vary by almost two orders of magnitude. 

A broad (outflow) component can be clearly identified in the 
spectrum, but its hyperfine components are not resolved. Thus, 
it is not possible to determine an excitation temperature of the 
broad emission from the spectrum. Assuming Tex - 100 K, a 
column density of A^oh = 2.0 x lO''* cm"^ is obtained by fitting 
the line wings only. A lower limit on the OH column density of 
A^OH ^ 7.2 x lO'-' cm"^ is derived by assuming that the excitation 
temperature is equal to the upper level energy (270 K) of the 
transition. Both fits yield a line width of Ao = 21.9 km s ' and 
a line position wisr - -35.7 km s"^ The line optical depth is 
T < 0.01 for these excitation temperatures. An optical depth of 
T = 1 is reached for Tex < 34 K and no good fit can be obtained 
for Tex ^ 32 K. In the following, Tex - 270 K is assumed. 

To convert the column density into an OH abundance in the 
outflow, the '^CO(3-2) observation of Hasegawa et al. (1994) 
was used to derive A^ ., = 1.3 x 10^' cm"^ with RADEX 
dvan der Tak et al.ll2007l) . assuming a CO/H2 abundance ratio of 
10-"*, a density of 10^ cm ^, and a temperature of 60 K, as de- 
rived in their paper This H2 column density converts the lower 
limit of A^oH ^ 7.2 x 10'^ cm"^ into a lower limit on the OH 
abundance of xqh = A^oh/A^Hj > 5.5 x 10"'^. 

3.2. Spherical model of the envelope 

In addition to the slab modeling, the narrow emission component 
is compared to the result s from full radiative transfer mod els us- 
ing the 'RATRAN' code (iHogerheiide & v an der Tak'2000). The 
physical structure of W3 IRS 5 is taken from van der Tak et al. 
(1200 0) with a power-law density profile ranging from « = 10^ - 
10^ cm"^ and temperatures 50 < T < 950 K within the Herschel 
beam at 1837 GHz. Dust and gas temperatures are assumed to be 
equal. Far-infrared excitation by the dust continuum for Tdust is 
included. We also adopt their distance (2.2 kpc) for consistency. 
Dust opacities based on grain s with thin ice mantles are assumed 
dOssenkopf & Henningll994i table 1, column 5). The molecular 
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Fig. 2. Resulting HIFI spectrum after subtraction of the best 
fit outflow component. Overplotted are the spherical envelope 
models for constant OH abundances of 1 x 10 (blue dashed), 
8 X 10"^ (red solid with individual components red dotted) and 
5 X 10 (pink dashed). Hyperfine components were simply 
added (see text). 

abundance is assumed to be constant throughout the envelope, 
because this is the simplest structure yielding a good fit to the 
line profiles. Figure|2]shows the comparison between model and 
data, where the best fitting slab model for the broad component 
(cf. Sec. 13. II ) is subtracted from the observation, leaving only the 
narrow component for a simpler comparison. 

Models assuming an OH abundance of xqh - (0.5 - 1) x 10"^ 
agree reasonably well with the observed narrow components. 
RATRAN does not treat line overlap and overlap effects can 
therefore not be treated accurately for optically thick lines. Thus, 
we have simply added the intensities of the components as would 
be appropriate in the optically thin limit. 

The spherical envelope models are also consistent with the 
79, 84 and 119 jum lines being in absorption, as observed in the 
unpublished spectral scan obtained with PACS on Herschel. The 
79 and 119 jim transitions are connected to the ground ^03/2 
level and therefore prone to absorption. A quantitative analysis 
is not possible because the central spatial pixel of the detector is 
saturated. 



4. Discussion 

Insight into the water chemistry can be gained by comparing OH 
and H2O abundances, because thes e species are linked th rough 
the OH -H H2 « H2O + H reactions. lChavarria et"aD(l2010l) iden- 
tify broad and narrow components in H2O fine profiles from 
W3 IRS 5 similar to OH. OH/H2O abundance ratios can be de- 
rived separately for gas of the envelope and the outflow. 

From the spherical non-LTE envelope models, which do not 
consider overlap of hyperfine components, we find an OH abun- 
dance of JCoH ~ 8 X 10"^ for W3 IRS 5. New calculations of 
the H2O abundance xh,o in the envelope, based on the work by 
Chayarn'a et al. (2010) but recalculated using the same physi- 
cal structure as for OH, yield x^^o - 10"^ for T > 100 K and 
■^HtO = 10"^ for T < 100 K. This allows us to estimate the 
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OH/H2O abundance ratios: H2O is about three orders of magni- 
tude more abundant in the inner envelope (T > 100 K), while the 
OH/H2O ratio is around unity in the outer part (T < 100 K). The 
same conclusion is reached i n the best fit H2O abundance struc- 
ture of lBoonman et alj (120031) to ISO and S WAS data, which was 
derived with the same physical model as adopted here. 

The inferred OH/H2O ratios in the envelope are consistent 
with the current picture of the water ch emistry in dense clouds 
(for detailed discussion o f processes see iHoUenbach et ani2009t: 
Ivan Dishoeck et aDl20I Ih . Different paths can increase the gas- 
phase H2O abundance in the inner part of protostellar envelopes. 
At r > 100 K, H2O starts to evaporate from the ice mantles of 
dust grains. When temperatures of T > 230 K are reached, which 
are typical for the innermost parts of high-mass protostellar en- 
velopes, H2O can be rapidly formed in the gas phase through 
OH -H H2 => H2O + H. In this case, OH is only a transient step 
of the H2O formation process and thus OH/H2O « 1 in the ab- 
sence of UV photons. In the outer envelope , low-temperature 
gas-phase chemistry provides jcqh ~ 10'^ (iDotv et alJ I2002L 
Fig. 1 1). At the outer edge, H2O and OH can be released into the 
gas p hase b y photode sorption from grain ma ntles. Obe rg et aTl 
(I2OO9I) and lAndersson & van Dishoeckl (l2008l) find from labora- 
tory data and theoretical calculations that roughly equal amounts 
of OH and H2O are released at low temperatures. Depending on 
the optical depth of the lines, this outer layer may dominate the 
emission. 

The slab model method outlined in Sec. 13.11 is also used 
to estimate the H2O abundance in the outfl ow component. 
Availability of HIFI para-H20 2o,2 - li,i (from IChavarria et al.l 
I2OIO ) and unpublished 2i 1 - 2o,2 (Chavarria et al. in prep.) and 
33,1 - 4o,4 data allow us to constrain A^p-HjO and the para-H20 
excitation temperature in the outflow simultaneously. We find 
A^p-H,o = 6.4 X 10"* cm-'- and thus A^h.o = 2.6 x 10'^ cm^^ 
for an assumed ortho- to para-H20 ratio of 3:1. The water abun- 
dance is calculated to be 2.1 x lO"*" with A^H2 = 1.25 x 10^' cm"^. 
Under the assumption that the OH and H2O emission arise from 
the same gas, we can calculate a lower limit of 0.028 on the 
OH/H2O abundance ratio in the outflow of W3 IRS 5. This 
limit is consistent with the emission origi nating in either a fast 
(v > 60 km s dissociative J-type shock (Ne ufeld & Dalgarna 
|1989|) or a slowe r UV-irradiated C-type shock (v ~20-30 km 
s~'. IWardlell999l) . Standard C-type shocks underproduce the ob- 
served column density ratio by one order of magnitude or more 
(Kristensen et al. in prep.). It is not possible to distinguish be- 
tween the two types of shocks with the current data and avail- 
ability of model results. 

5. Conclusions 

The OH fines at 1837.8 GHz from W3 IRS 5 detected wifli 
HIFI consist of two components: a broad component, where the 
triplet components are blended because of the large individual 
line widths, and a narrow component with velocity resolved hy- 
perfine structure on top. These results indicate that OH emis- 
sion observed at a lower spectral resolution, e.g. with PACS, 
can be a blend of envelope and outflow contributions. Low- 
mass sources reach lower temperatures than the high-mass coun- 
terparts at the projected distance of the Herschel beam. From 
RATRAN models of low-mass sources we thus find a lower ex- 
citation, resulting in a weaker envelope contribution to the OH 
1837 GHz lines than from W3 IRS 5, consistent with the results 
for HH 46 in Wampfler et al. (2010). Because the broad com- 
ponent also appears in the line profiles of other species from 
W3 IRS 5, in particular those of CO and H2O, the outflow is 



the most likely origin of the emission. Comparison with H2O 
yields a lower limit on the OH/H2O abundance ratio of 0.028, 
consistent with an origin from a fast J-type shock or a slower 
UV-irradiated C-type shock. The narrow emission agrees well 
with spherical envelope models with a constant OH abundance 
of xqh ~ 8 X 10"^. Co mparison with new H2O results based on 
Chavar ria et aP (1201 Ol) give OH/H2O 1 for T < 100 K and 
OH/H2O 10'^ at r > 100 K, consistent with the current pic- 
ture of the dense cloud chemistry with freeze-out. 
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Appendix A: Gaussian fitting 



Table A.l. Gaussian fit results to the OH line components using 
the velocity range [-100,30] and first order baselines. The v - v\st 
velocity scale is given relative to the laboratory frequency of the 
2+ 1- component = -38.4 km s^' = 1837.817 GHz). 
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Fig.B.l. Comparison of the broad component of OH (black), 
H2O 2o2 - 111 (red), and CO(10-9) (blue) observed with HIFI. 
The black dashed line indicates the source velocity (uisi - 
-38.4 km s-'). 



B of iBruderer et al.l (l2010l) . Both OH line components (nar- 
row/envelope and broad/outflow) are represented by a slab in 
front of a continuum source. The continuum temperature Tcom is 
derived from the observed single-sideband continuum value and 
the source is assumed to be fully covered by both slabs. No ge- 
ometry is included, except that the outflow slab is placed in front 
of the envelope slab. Each slab has four free parameters per line: 
the OH column density A^oh, the excitation temperature Tex, the 
line width Ad, and the position (in velocity space) of the line cen- 
ter Disi-, but the excitation temperature is assumed to be the same 
for all hyperfine transitions. The normalized level populations of 
the upper (Xy) and lower level {x\) are therefore determined by 
the Boltzmann distribution at Tex, 



" Parameter fixed. The given errors are only fitting errors and do not 
include the calibration uncertainty. 



Appendix B: Line width comparison 

The comparison of OH and H2O column densities in the out- 
flow is based on the assumption that the emission arises from the 
same gas. Figure IBTI illustrated the similar widths of the broad 
components of OH (Av ^ 27 km s '), CO {Av 29 km s">), 
and H2O (Ay 26 km s '). The broad component of OH is a 
blend of three hyperfine components with individual widths of 
Av ^ 22 km s"', derived from the best fit slab model (Sec. 13. U . 



— - — exp 

Xl gl \ ksTe^ 



(C.l) 



with the statistical weights and g\ of the upper and lower 
level, respectively, Boltzmann's constant ^b, Planck's constant 
h, and the frequency vq of the transition. The radiation temper- 
ature T (v) is derived from the solution of the radiative transfer 
equation, using the Rayleigh-Jeans approximation, as 



T{v) = reo„t(v)exp(-T'=-)exp(-r°«) 

+ 7^Bv„(r|r)[l-exp(-T"^)] 

with c being the speed of light, By^ the Planck function, r''™ the 
optical depth of the envelope layer, and the optical depth of 
the outflow layer The line optical depth of every slab is the sum 
of the contributions of all M hyperfine components and can be 
calculated from 

*f 2 / \ 
-^^)=^oh|:^A^./|-.^),^(v) (C.3) 



where Aui is the Einstein coefficient and (p(v) the normalized 
line profile function of the transition, assumed to be a Gaussian 
of width An centered at uisr. This approach takes line overlap into 
account. 



Appendix C: Slab modeling details 

The modeling of the OH line spectrum (cf. Sec. 13. Il l is car- 
ried out with the slab model code presented in appendix 
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Fig. C.l. Icr and 2cr contours for the fit of the narrow (envelope) 
component with slab models. The best fit is indicated by the red 
dot. 



